Novel Geometrical Models of Relativistic Stars. 
II. Incompressible Stars and Heavy Black Dwarfs 
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In a series of articles we describe a novel class of geometrical models of relativistic stars. Our 
approach to the static spherically symmetric solutions of Einstein equations is based on a careful 
physical analysis of radial gauge conditions. 

It turns out that there exist heavy black dwarfs: relativistic stars with arbitrary large mass, 
which are to have arbitrary small radius and arbitrary small luminosity. In the present article 
we mathematically prove this new phenomena, using a detailed consideration of incompressible 
GR stars. We study the whole two parameter family of solutions of extended TOV equations 
for incompressible stars. This example is used to illustrate most of the basic features of the new 
geometrical models of relativistic stars. Comparison with newest observational data is discussed. 
PACS number(s): 04.20. Cv, 04.20. Jb, 04.20.Dw 



I. INTRODUCTION 

This is the second of series articles in which we de- 
scribe a new geometrical models of general relativistic 
stars (GRS). One can find the general scheme, notations, 
basic equations, additional conditions and basic princi- 
ples and properties of these models in the first article of 
this series - [1]. The essentially new features of GRS in 
these new models are not based on the critics or revision 
of very general relativity (GR). They are a result of more 
deep understanding of its applications, and on solution 
of some open problems in this theory, like the physical 
justification of the choice of GR gauges. The preliminary 
knowledge of the article [1] is highly recommended. It is 
essential for the right understanding of the present one. 

Here we consider the simplest specific example of ap- 
plication of the general scheme, described in [1]: GRS, 
made of an incompressible matter with equation of state 
(EOS): 



e{p) = const 



: e. 



(I.l) 



In Eq. (1. 1 ) £ is the energy density and p is the pressure of 
stelar matter. We are using units c = Gat = 1- Further 
on the letter "C", as an index of different quantities, 
denotes their values at the center of the star, the sign 
" * " - the values at the edge of the star. 

Such simple model has a limited physical significance, 
because it leads to an infinite speed of the sound in the 
fluid [2]. Nevertheless, its consideration is useful, be- 
cause: 

• It gives some general estimations for the properties 
of the solutions of extended Tolman-Openhcimer-Volkov 
(ETOV) equations with arbitrary EOS [2]. 
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• This simple model can be used as a good physical ap- 
proximation for description of neutron stars with matter 
density w 2.85 x 10^^ g/cm^ and ^ 5 x 10'^^ erg/cm^, 
because under these conditions the nuclear matter be- 
haves much like incompressible fluid [2]. 

• The EOS (I.l) has the advantage that it yields an 
analytically solvable model of stars. The simplest degen- 
erate set of solutions with luminosity variable pc — was 
at first found in the Schwarzschild pioneer article [3] . 

• It is a basic example of application of GR to the 
stelar physics [2]. 

Because of the overestimation of the role of the lumi- 
nosity variable p, up to now the considerations of GRS 
structure were based, as a rule, on the Hilbert gauge 
(HG): p(r) = r, in which the condition pc = seems to 
be natural. Here r is the proper radial variable of the 
spherically symmetric problem at hand [1]. 

To the best of our knowledge, the solutions with pc > 
have not been studied and do not have a proper phys- 
ical interpretation. We intend to fill this gap in the 
present article. The general solution of this problem, 
described here, turns to be much more complicated and 
more interesting then the degenerate one, considered by 
Schwarzschild in [3]. 

The new models of GRS recover an essentially new and 
unexpected relativistic physics. In particular, in these 
models GRS with arbitrary large mass m* are allowed. 
These are to have arbitrary small geometrical radius i?* 
and arbitrary small luminosity. We refer to such amazing 
relativistic objects as heavy black dwarfs (HBD). 

In the present article we give a mathematical proves 
of the existence of HBD in the specific case of incom- 
pressible matter. Taking into account that the EOS has 
a quite week influence on the proper radial gauge p(r) 
[1], one may expect this result to be general. 

Pure GR reasons, which are independent of EOS, can 
not yield restrictions on the maximal mass of stars. Con- 
straints of that kind may arise only due to quantum 
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statistics via the EOS, and in equal footing both in New- 
ton gravity [2, 4] and in GR. These theoretical conclu- 
sions may give a more precise physical understanding of 
the real phenomena in stelar physics and need a further 
study. 

The new family of GRS is two parameter one. As a 
free parameters one can consider, for example, the mass 
m, and the radius FL, [1]. Between these quantities, in 
general, there exist no stiff functional relation, as in the 
widely accepted GR models with extra condition pc = 
[2, 5]. Actually, just this assumption is responsible for 
existence of the well known mass-radii relations [1]. In 
a framework of a given theory of gravity (Newtonian, 
GR, etc) their specific form depends on EOS, but their 
common origin is in the condition pc = 0. 

The very recent observational data (see J. Madej et 
al., 2004 in [5]) seem to confirm the existence of two- 
parameter family of white dwarfs with free parameters 
m, and in proper domain, thus rejecting the existence 
of stiff mass-radii relations in Nature. Here we give a pre- 
liminary comparison of these data with our new models 
of GRS, using the results, obtained for the simple case of 
incompressible stars. 

Making use of the simplest EOS (I.l) we will be able to 
illustrate in most transparent way the differences between 
our approach to the GRS structure and the commonly 
accepted one with extra condition pc = 0. 



and 



dip m + inp^p 



>0, 



dp p{p — 2m) 
v{pc,pc) = ^c, 'p{p*,pc) = 'P* 

have the form 

p 

wo(p,Pc) = 47re / p^^J-gpp{p, Pc)dp, 



(11.6) 



(11.7) 



^(P,pc)=ln l±(^?£iM£l^U^c. (II. 



\/-9pp{p, Pc) 



In the above formulas 



P 3 

x(P,Pc) = 47rir j p(^\J -gpp{p, pc)^ dp. (II.9) 



PC 



The solution (II. 4) can be rewritten in a more convenient 
form 



X(P, Pc) = 



V 9pp ( V 5pp 



w \ w + 1 



(11.10) 



II. GENERAL SOLUTION OF ETOV 
EQUATIONS FOR INCOMPRESSIBLE MATTER 
IN HILBERT GAUGE 



For EOS (I.l) the equation 



dm 
dp 



A-KEp^ > 0, 



(11.1) 



miPcPc) = 0, m{p*,pc)=m^ 



splits from the whole ETOV system of differential equa- 
tions and has a simple general solution: 



m{p,pc) = gTTE (p^-Pc) • 



The general solution of the equation 

dp {p -|- e) (to -|- 4:Trp'^p) 
dp 



<0, 



p{p - 2m) 
p{pc, Pc) = PC, p{p*,Pc) = 
can be written in the form 



(II.2) 



(II.3) 



/ X , \/-9pp{p,Pc) {pc+s) „j .s 

p{p, pc)+£= —-, , ^-.s , r ■ II.4 

1 + {Pc/£ + l)x(/0,Pc) 

The general solutions of the equations 



dmo 

dp 



= ^T^sp^^f-app > 0, 



(II.5) 



moipcPc) = 0, mo{p*,pc) = mo*, 



making use of the scale-invariant ratio w = p /£[!]■ 
Then the metric coefficients can be written in the form: 

gu{p,pc)=e'^^''''"''\ 

gpp{p,Pc) = z — ^T—f ri-. (11-11) 

l-2m{p,pc)/p 



III. SCALE PROPERTIES OF THE PROBLEM 
AND HG SCALE INVARIANT QUANTITIES 
FOR INCOMPRESSIBLE STARS 

The EOS (I.l) is not scale- invariant. Since the scale 
factor is A = const, the scale transformations: 

p^Xp, m^Xm, too^Atoq, s^X~^e, p^X~^p. (III.l) 

described in [1, 6] in more details, map the solution of 
ETOV system with a fixed value of the constant density e 
onto another solution of the same system with a constant 
energy density e' ~ X^^s. Utilizing this property one can 
consider only the case e = 1. Instead, we prefer to use 
a scale invariant variables. In contrast to other authors 
[6] , we avoid the a prior choice of such variables and will 
find them during the very course of the solution of ETOV 
system with EOS (I.l). 

In the case of incompressible star one can use the rep- 
resentation 



9pp{p,Pc) = -^P^ /Pi{p,Pc), 



(III.2) 
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where 



and 



is a naturally normalized polynomial of fourth degree, 
which roots are /3i,2.3,4- It is easy to find these roots in 
the form 



1 



Pi 



P-i 



^2TTe (1 + cr2/3) 
^27re(l + cr2/3) 



> P2 = > 



> Pi 



4(1 + ^) 



V27re (1 + ct2/3) ■ 



(III.4) 



where the A-invariant parameter cr is defined via the sub- 
stitution PC - ^^^^^^^^.^3^ (I^)^^^IO]. 

Now it becomes clear that instead of the luminosity 
variable p it is better to use the A-invariant one, ^: 



e V2vre (1 + a^/S) p = W 1 - p, (III.5) 



Stte 



where G (0,1) is the only positive root of the cubic 



equation 



3 c2 



As a result we obtain 



3 ^ ? 



(III.6) 



The naturally re-normalized A-invariant polynomial 

4 

^c) =n (^-^^) =^(e - 1) (fH+ec) (ni.7) 



ei-i>e2=o> 



^3 — 7i 



>e4 



if ec e [1/4, 1] 



(IIL9) 



There exist two degenerate cases: 

a) Schwarzschild degenerate case: = = ^3 = 
is a double root and P4(C,0) = -^^(1 - C^); 

b) A new degenerate case: = \/l/4 =^ ^3 = = — ^ 
is a double root and ^4(^,0) = -^(1 - 0(C - 1/2)^ 

In both cases the corresponding elliptic integrals, 
needed for an explicit solution of the problem, are re- 
duced to elementary functions - see Appendix A. 

For the A-invariant quantities 



Po{£.,£.c) = 



Stte 



mo{p,pc) (III.IO) 



from Eq.(II.2), (II. 7) we obtain (see Appendix A) 



Mo(C,^c) 



21-er 



V 



^={j3ll{^,^c) - J3|l(Cc,?c)). (in.ll) 



has the following re-scaled roots: 

Polynom 





\ \ \ 








\ ^-^^\\°-'' 
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FIG. 1: The polynomial P4,{^,^c) for different values of the 
parameter ^c- 



6 = 1 > 6 = > 

?3 > 44 - 



i+yr^i^ 




FIG. 2: The function p{^;^c) £ [0, 1/2] and the lines ^ = 
const and = const. 

Instead of the local binding energy Am{p, pc) '■= mo — 
m, which is not A -invariant, one can consider the ratio 



if [0,1/4], 



(III., 



g(p, pc) := m/mo = m/Mo = git Cc) e (0, 1). (III.12) 



c 



FIG. 3: The function /io(Ci £,c) and the hnes ^ = const and 
5c = const. 



FIG. 4: The function x(Ci Cc) and the hnes ^ = const and 
5c = const. In contrast to mo(5;5c'), the function x(C)Cc) is 
not bounded from above. 



It measures in a A -invariant way the local mass defect 
of the star mater, i.e. the mass defect in the sphere with 
luminosity radius p and center C. 

Another important A-invariant local (in the above 
sense) quantity is f{p,pc) = q{p, PcY ~ 9tt{p, Pc) + 1- 
In the case at hand it has the form 

f{p,Pc)^(—)\ — = 9' + e, (III.13) 
\moJ p 

where = ^ > is the local compactness of the star. 

For incompressible stars — ^1^%^. 

Using the results (A. 5) and (A. 6) of Appendix A, after 
some algebra one obtains the expressions: 



_c_ 

I J y^iWM 

t,c 



'e-ec 



{Mi{^,^c)-Mi{^c,^c)) (111.14) 



and 



V^l-e^(^4|3(C,?c) - J4|3(ec,&)) (111.15) 



3 
2 

for the basic A-invariant local quantities g and x- 

As one sees, the technical problem of finding the gen- 
eral solution of ETOV system for incompressible matter 
in the most general case is reduced to the calculation of 
the elliptic integrals in Eqs. (III. 10), (III. 14) and (III. 15). 
This calculation is described in the Appendix A. 



IV. RADIUS OF INCOMPRESSIBLE STARS 

Now we are able to proof the existence of a finite 
coordinate radius < oo of incompressible stars. It 
corresponds to some finite value of the luminosity vari- 
able ^* < oo and to some finite geometrical radius 
i?* = 3l*/v/87re(l-C3,)/3 < oo. (Here 3^* is the di- 
mensionless A-invariant geometrical radius of the star.) 

By definition, in A-invariant terms the edge of the star 
is defined as a point ^* at which 



p(^*,Cc)=0. 



(IV.l) 



Proposition 1: For all solutions of ETOV equations 
(II.2)-(II.6) with EOS (I.l) and arbitrary £,c > there 
exist a unique solution £ (^Ci ?*'^**(Cc)) of the Eq. 
(IV.l). Here the value C*"*(Cc) € (^C) 1) corresponds to 
the nonphysical limiting solution with infinite value of the 
pressure at the stelar center: p^f** = oo. 

Proof: 

1) Existence of unique solution of Eq. (IV.l): 
From Eq. (III. 6) and (III. 15) one easily obtains in the 
limit ^ 1 - 0: 



{2+ecfi-i 



Hence 



lim 



^^^^^^.^e (2/3,1) 



(IV.2) 



(IV.3) 



5 



As a result 

Mm wi^,^c)^-^—^<0, (IV.4) 

i.e. the continuous, strictly monotonic function w{£,,£,c) 
decreases on the interval [^C; 1] from some value wq = 

w(Cc,Cc) > to the value w(l,Cc) = < (i-e-, 

w{^,^c)\[wc-{l-ec)/i] for e G Kc,l)). Then there 
exist a unique value S (^c", 1), such that , Cc) = 0. 
This value ^* defines the radius of the star. 

2) Limitations on the A-invariant radius of star 
The Eq. (11.10) permits us to rewrite the definition of 
the stelar radius (IV. 1) in the following explicit form: 




FIG. 5: The function wc{£,*',^c) and the lines ^* = const 
and = const. Note that wc(C*;Cc) = oo on some special 
line 5* = 5t'^**(^c) which will be studied in details in Subsec- 
tion E of the present Section. 

It demonstrates the relation between the value wc at 
the center C of the star and its radius This relation 
is the main difference between relativistic theory of stars 
and Newtonian one. 

Using the following properties of the functions 
\/^9pp{^*,£.c) and x{^*,^c) 

a) - 9pp{£.c,^c) = 1; 

b) V-9ppii*^^c) / [l,oo), x{i*.ic) / [0,oo) - for 

c) forO<e-Cc^2^y5We^-x(e*,?c) >0; 

d) fore^l: V-5pp(C*,ec)/x(^*,Cc)«(2 + ^^)/3e 
(2/3, 1); , one easily obtains that: 

i) wc^icic) = 0; 

ii) wc{l,^c) = -1; 



ii) There exist a unique point ^J''** € {^c, l)j such that 
v/-5pp(e"*,Cc) = X(e"*,ec) => wciiT' - 0,Cc) = 
+00 and wc(Cr'* + 0, ic) = -oo. 

This means that there exist a limiting nonphysical 
solution - an incompressible star with a critical radius 
^^"'^ G {£,c, 1), which corresponds to an infinite pressure 
pg"'* = oo at the center C. Hence, the A-invariant radius 
of the incompressible stars is constraint in the interval 
{^c,£,Z^*)- It is obvious that is a function of S^c- 
At the end of this Section we shall obtain the equation 
for determining of the function £,^^'^*{£,c) and study its 
solution in details. 

This completes the proof of our Proposition. 

The luminosity radius of the star can be find in the 
form /9* = l^pc > PC and is not bounded from above if 
e 0, or/and for ^ 1. 

Having in our disposal the A-invariant radius i^* of the 
star, we are able to introduce another basic A-invariant 
characteristics: 

i) The total Keplerian mass = m(C*7Cc) G (0,1/2) 
- for < < < 1, see Fig 2. 

ii) The total proper mass /io* — ^J■a{£,*y^c) > 0. 

iii) The total mass defect ratio — q{£,*,^c) > 0. 

iv) The total compactness = <;(C*;Cc) = 2^1^ /p^ = 
j^i^e(0,l), etc. 

v) In addition we obtain the relations 

ip,=ip{^,,^c) = l^iwc + l) + ipc (IV.6) 
and (see Appendix A): 

yi^(ji|i(e,Cc)- Ji|i(^c,Cc)). (IV.7) 

It is easy to see that 3?*(^*,^c) — > oo when — and 
^ 1 - 0. 

As seen in Fig. 2 the function p,^,{S,*,£,c) has values in 
the interval [0, 1/2]. At the same time we obtain a new 

Proposition 2: In the specific limit: — > 1, — ^ 1, 
< !i ^* ^ 1 we have 

lim /i,(C*,Cc) =:p*(l-0, 1-0) e [0,1/2], (IV.8) 

5c<5.<i 

i.e., the limit /i*(l — 0, 1 — 0) is bounded, but not definite 
and can have any value in the interval [0, 1/2]. 

The analytical proof becomes obvious from the rep- 

resentation ^*(^*,^c) = 2 i+gc+C consider- 

ation of the limit ^c" ^ 1 of p.*{£,*,£,c) on the curves 
i* = £,c + k{l — ^c) with an arbitrary k S [0, 1]. 
Note that for any k € [0, 1] and rt > l,p > 1: 

i) On the curves = -|- k{\ ~ the limit (IV.8) 
has an universal value 0. 

ii) On the curves = 1 - A:(l - C^)" the limit (IV.8) 
has an universal value 1/2. 



6 



R 




FIG. 6: A part of the function Ji, ; ^c) and the hnes ^, 
const and = const. 



Hence, the typical behavior of the hmit (IV. 8) is given 
by the last two cases: i) and ii). The case of linear de- 
pendence between and is an exceptional one. 

The simple, but important property of the function 
fJ'*{^*,S.c)i described in Proposition 2, will influence fur- 
ther results in the theory of incompressible relativistic 
stars, when the same limit of other quantities will emerge. 



V. THE BASIC Rm MAPPING FOR 
INCOMPRESSIBLE STARS 



The basic mapping 



Rrn 



U*(3?*,M*),^c(3?*,M*)}. (V.l) 



was defined in general case in [1]. Here we study this 
mapping in the specific case of incompressible GRS in 
scale-invariant variables: 

Substituting in the formula (IV. 7) ~ f*(M*iCc), or 
= ^cip*T^*)y obtained from the first of the relations 
(III. 11) in the form 



(p* , Cc) = ^/2iI7+{l^^2iI^, (V.2a) 

Cc(/^*,C*) = VW^^^i^^^W^^^ilT), (v.2b) 

one obtains the functions 3?* = ^*{fJ-*,S.c) and 31* — 
X(Ai*,C*)- 

The corresponding inverse functions = (/i* , 31^ ) 
and = Cc(/^*, 3?*) define the basic mapping Rm (V.l) 
in A- invariant terms. This mapping is illustrated in fig- 
ures 8 and 9. 

As a byproduct from Eq. (V.2b), 2/x* < 1 and > 
we obtain a new inequality > 2/i*. Then, since < 




FIG. 7: The function = ?c(m*,^*) (V.2b) for different 
fixed values of parameter G (0, 1/2). 




FIG. 8; The function ^.(^.;K,) G [0,1]. The behavior of 
this function in the limit /i, 0, 3?* — > is similar to the 
one, described in Proposition 2 for the function /i*(^, ,^c)■ 



^* < 1, we have 



2u* 2u* 



(V.3) 



For incompressible stars the relations, which determine 
the domain d[^^^^ [1], acquire the following specific form 
(For the used notations see Appendix A.): 



< < < 1, 

-e*(i + e*) < o< ^i^. 



(V.4a) 
(V.4b) 



1 



3 



< 



2(^4|3(e*,ec))- J4|3(ec,^c))< 



(V..4c) 
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As a result of the first of the relations (V.4) the con- 
dition (V.4b) is fulfilled. 

(2) 

Parts of the corresponding domains are shown 

in Figs. 10 and 11, where the condition (V.4c) is still not 
taken into account. 



D.5 




FIG. 10: The domain D^^J^^ for the function ^c{fi*;^*)- 



VI. THE CRITICAL RADIUS AND CRITICAL 
MASS OF INCOMPRESSIBLE STARS 

The limiting case of Eq. (V.4c) yields the following 
equation for the critical value S,^"*^ = £,!^"'^{^c) of the 



dimensionless luminosity variable: 

o tcrit 

.«{e-,&)-./.,(fc.fc)^5^^=^=. ,v,.i, 

It is obvious that this equation defines a new function 
^*"'^{£,c), which describes the basic difference between 
GR and Newtonian incompressible stars. There are no 
any physical parameters in the Eq. (VI. 1). Hence, the 
function ^^"'^{^c) is an universal mathematical one. 

In the degenerated Schwarzschild case when £,c = 
we have — -P4(C,Cc) = — ■f^) and the integral in Eq. 
(VI.l) can be easily done. This gives (see Appendix A) 

1 _ 2 

Hence, 

C;"*(?c = 0) = ^/8/9 « .942809. (VI.2) 

Thus we have obtained the exact value of the new func- 
tion C*'^'*(Cc) at the point — 0. It is clear that 

= 1) = 1. 

The shape of function S,^"*{^c) is shown in Fig. 12. 
Some of its basic values are given in Appendix B. As 
seen, the values of the function CJ'^**(Cc) are in the inter- 
val [a/8/9, 1]. Hence, for a given value of A-invariant lu- 
minosity variable > 0, the value of ^* varies in the in- 
terval Kc, CS"*(ec)], with upper limit ^r\^c) > yW^- 

This changes radically our understanding of the rel- 
ativistic theory of stars, because now in the domain 
< < 6 < 1: 
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0.2 0.4 0.6 



FIG. 12: The function 5r"(Cc) as a border of the physical 
domain of the variable f« € (^c^r^C^c))- 



6 
5 

4 

m.3 



m 



crit 



Physical domain 
02 0.4 ^ 0.6 08 1 

FIG. 13; The critical mass of the star m!;;'^*(Cc) as a border 
of the physical domain of the variable m* € (p,mf^*{^c))- 



As a result of the last inequality one obtains the famil- 
iar relativistic restriction for the Schwarzschild degener- 
ate case: 



(2m,/p* ) 



Uc=o 



This constraint yields the following restrictions on the 
luminosity variable and Keplerian mass of Schwarzschild 
incompressible stars: 



1 4 1 

m*| < - . 

'Sits \sc=o " v37re 



(VI.4) 



Hence, in the relativistic model with pc = we are not 
able to describe stars with fixed density e = const and 
with arbitrary large mass . 

In contrast, in our general geometrical model of in- 
compressible relativistic stars with arbitrary fixed value 
of luminosity variable € [0, 1] we have the restrictions: 



1 



e"*(ec) 



00 : for ^c— > 1 — 0, 



oo: (VI.5) 



for ^ 1 - 0. 

These upper limits pS'''*(fc) and mj'"'*(^c) are ea;act, i.e., 
the corresponding quantities p* and m* can become ar- 
bitrary close to them for a given value of ^c- 

The calculation of the limit of p';"*{^c) is obvious from 
^^"\^c = 1-0) = 1, see Fig. 12. The analytical proof of 
the above limit for m^''**(^c) when — * 1 — follows im- 
mediately from our Proposition 3 (see Section IV, C) and 
Eq.(III.lO). The behavior of the critical mass ml"*{^c) 
as a function of the central value of luminosity vari- 
able is shown in Fig. 13 for density 327r£/3 = 1 [11]. 

Thus we have arrived at the following 



Proposition 3: For any fixed density e = const in 
our essentially non-Euclidean general model of relativis- 
tic stars there exist incompressible stars with arbitrary 
large mass to* . 

Obviously, stars with arbitrary large mass m* exist 
only for large enough values of luminosity variable ^c- 

From Eq. (IV. 7) one obtains the following critical 
value for the geometrical radius of our relativistic incom- 
pressible stars: 



Jiii(er*(^c),ec)- Ji|i(gc,$c 

i/87r£/3 



(VI.6) 



The behavior of the critical geometrical radius RZ^*{^c) 
as a function of the central value (^c of luminosity variable 
is shown in Fig. 14 for density 327r£/3 = 1. 



R., 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 



Physical domain 



cnr 



^] 0.2 0.4 0.6 0.8 1 

^c 

FIG. 14: The critical geometrical radius of the star il*'^**(^c) 
as a border of the physical domain of the geometrical stelar 
radius € (il^"*(^c), oo). 

Thus, we have proved the existence of new amazing 
regular relativistic objects with constant matter density: 
They can have an arbitrary large mass to*. As a result 
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of strong gravity, their geometrical radius i?* and lumi- 
nosity ~ sre arbitrary small, when the mass 
is large enough. At the same time in presence of such 
objects the geometry of the space-time is regular every- 
where. Horizons of any type do not exist. 

It seems reasonable to refer to objects as to heavy black 
dwarfs. As seen in Fig. 13, the existence of HBD is 
impossible under widely accepted extra condition pc = 0. 



The critical proper mass 



'(^c) is bigger then 



i^T i^c)- Hence, it goes to infinity in the limit ^ cx3 
together with 'm';"*{^c)- 

More interesting characteristic is the A-invariant criti- 
cal mass defect ratio e^''**(Cc)j obtained from Eq. (HI. 14) 
in the form: 



(VI.7) 



3(J3|l(e"*(^c),Cc)-^3|l(ec,ec))^/T^ 
This function is represented in Fig. 15. 
1 



VP, 



0.8 
0.6 
0.4 
0.2 





Physical domain 



0.2 0.4 0.6 0.8 



FIG. 15: The critical mass ratio gT^^{£,c) as a border of the 
physical domain of the variable g, £ l) ■ (In the 

Fig 15 the symbol Vp stands for q.) 

As seen, the critical value of the mass ratio g';"'^{^c) 
decreases monotonically from 



Ifi / /Q 

er*(0) = y 1-3^-arcsin 




.609477 



to 1/2 when the variable ^c" increases from to 1. 

Hence, in the case of incompressible mater the ex- 
tremely strong gravitational field, which arise in the limit 
i?* — > 0, m* — > c» is able to extract no more then one 
halve of the initial bare mass mp* of the stelar matter, al- 
though the relative part of gravitational energy increases 
in this limit, together with the mass defect. It can be 
shown that this limitation is EOS-dependent. 

One can understand the above behavior of the critical 
mass ratio Q^"*{(,c) analyzing the form of the functions 
f^r'i^c) - Fig. 16 and pf,:^\^c) " Fig. 17: 

A more deep understanding of the mass defect one can 
obtain looking on 3D surface of the function, £'(^*, S,c) = 



0.5 
0.4 
0.3 
0.2 
0.1 



,crit 



Physical domaiQ 



0.2 0.4 0.6 0.8 1 
^c 

FIG. 16: The critical function pT^^iS.c) as a boundary of the 
physical domain of the variable p, £ (^0, p'^^^^ {£,c)) ■ As seen, 
on the specific curve ^, = the function pT^^{£,c) 

goes to 1/2, when —* 1. (Compare this specific result with 
Proposition 2 and with comments after it.) 



1 




0.8 




,*0.6 




0.4 


Physical domain 


0.2 







0.2 0.4 0.6 0.8 1 



FIG. 17: The critical function /io"*(5c) as a boundary of 
the physical domain of the variable pot G (O, /io"*(Cc)) • 



fJ,//J, 




FIG. 18: The function g{^.,^c) = pi^„^c)/poi^„^c) > 
4/37r for incompressible star. The shadowed triangle is a part 
of the horizontal plane 4/3tt ^ .424413. 
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K^*,£.c)/MC*;^c) > 4/37r, shown in Fig. 18. The last 
Umit from below of the mass ratio g{C*,£,c) originates 
from EOS for incompressible stars Eq. (I.l) and will be 
different for relativistic stars with other EOS. It reflects 
the form of the right border of the 3D surface: 



£',(^*,0) = - 



3 arcsin — vT— 



(VI.8) 



calculated for the Schwarzschild degenerate case of in- 
compressible stars. An analogous simple expression can 

be drown for g y/ 1 /4^ , using results, given in Ap- 
pendix A. 

The real limit from below on f?(C*jCc) in the problem 
at hands is 1/2. It reflects the existence of the curve 
^i"*{^c), which is not shown in Fig. 18. 



R. 



1.2 
1.0 
0.8 
* 0.6 
0.4 
0.2 





Physical domain 



3 4 



FIG. 19: The function Rl"'' {mi ) as a boundary of physical 
domain for incompressible star with energy density ^tte = 1. 



VII. MASS-RADII RELATIONS FOR WHITE 
DWARFS AND THE NEW GEOMETRICAL 
MODELS OF GRS 

Our previous consideration [1] showed that there exist 
two different approaches to the GRS with stiff mass-radii 
functional dependence: 

1) The standard one, based on the extra condition 
p(0) = PC = 0, and 

2) A new one, with extra condition: F- fixed. It arises 
naturally in basic regular radial gauge [1]. Here F S 
(—00,00) is the stelar crust parameter. It defines the 
jump of the derivative of luminosity variable p{r) as a 
function of the radial one, r : p'(r* — 0) = e~'^p'(r* -|-0), 
at the edge of the star r*. 

At the same time one can consider a geometrical mod- 
els of GRS without any additional condition of that type. 
In this most general case stiff mass-radii relations, based 
only on GR considerations, i.e. independent on EOS, do 
not exist [1]. In this case m* and i?* are free parame- 
ters, constrained in some 2D domain, the precise form of 
which depends on the EOS. 

The physical domain of variables and i?* for in- 
compressible GRS is shown in Fig. 19. It is obvious that 
this physical domain is very similar to the one, recently 
observed for real white dwarfs, see Fig. 2 in the article 
by J. Madej, M. Nalezyty and L. G. Althaus in [5]. This 
demonstrates how we are able to reproduce the real ob- 
servational data in a qualitatively right way, using the 
simplest new geometrical model of incompressible GRS. 

The most important consequence of our analysis of the 
observational data is the conclusion that actually we do 
not observe stiff mass-radii relations in Nature. Instead, 
these data demonstrate a clear indication that m* and 
i?, are filling a 2D domain with a sharp boundary of the 
type, similar to the one, shown in our Fig. 19. This 
corresponds to our general models of GRS without any 
extra conditions [1]. Therefore we will skip the consid- 
eration of such relations in the present article, although 
they may be of some interest, too. 



For a more precise qualitative treatment of the obser- 
vational data one need to obtain new geometrical models 
of the GRS with realistic EOS. This may give a basis for 
explanation of the mass distributions of special type of 
stars, like white dwarfs, or neutron stars, independently 
of theory of gravity. We will present the corresponding 
results elsewhere. 



VIII. CONCLUDING REMARKS 

In conclusion several remarks have to be made: 

1. In the present article we have recovered an inter- 
esting new GR phenomena, studying the most general 
solutions of ETOV equations for incompressible matter. 
We were able to find proper physical meaning of all solu- 
tions of these equations, overcoming the requirement for 
regularity of solutions at the zero value of the luminosity 
variable p. As we have seen, this way we obtain a two 
parameters family of solutions. 

2. We wish to comment once more the most important 
result of present article: the proof of the existence of new 
GR objects - heavy black dwarfs (HBD). 

From physical point of view this becomes possible, be- 
cause the fulfilment of the condition p > pc > PG [1]) 
where pa = 2to* is the Schwarzschild radius, is equiva- 
lent to the introduction of non-permeable potential wall 
in the space of the luminosity variable p. Because of the 
Gauss theorem, this introduces in the energy balance of 
the body an effective "potential energy" with an infinite 
jump. 

Then, following the original Landau energetic consid- 
eration of the problem [9] , we conclude that the presence 
of any additional energy terms, due to EOS, can not de- 
stroy the existence of static equilibrium of the body, at 
least in the case of non-ultra-relativistic matter. This 
"p-space picture" explains in a more usual language the 
consequences of the new boundary conditions at the real 
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center of the star, under which we are solving Einstein 
equations. See for details [1]. 

In terms of the proper radial variable r no wall exist 
at all. Thus we see once more, that one is to give up the 
idea to consider the luminosity variable p as a measure 
of distances. It is responsible only for determination of 
geometrical area and luminosity of physical objects. 

3. The results of the present article raise many new 
questions and problems. 

For example, we obtain a new real alternative to the 
black holes (BH) for explanation of nature of the massive 
compact dark objects, observed in astrophysics. Hence, 
we need a new criteria to distinguish BH and HBD. 

An obvious problem is the possible relation of HBD 
with the observed could dark matter in the universe, as 
well as the relation of HBD with gravitational collapse of 
usual bodies. It is clear, that the present series of articles 
calls for a novel consideration of the collapse, too. 

For a comparison with observational data, we obvi- 
ously need consideration of new geometrical models of 
GRS with more realistic EOS for stelar matter. 

We intend to address all these questions in the forth- 
coming articles in this series. 
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APPENDIX A: CALCULATION OF THE 
ELLIPTIC INTEGRALS 

Using standard incomplete elliptic integrals of first, 
second and third kind {F, E and H) [7], we can write 
down the integrals 



which arc needed for the relativistic incompressible star 
problem, in the following form: 
I. n = 0, m = 1: 



=9-F{iz,k), 



In the generate cases one obtains: 



(A.2) 



a) Jo|i(e,0) = iln 

b) Jo|i(^,i/l74) = ^arctan(^^Il^ 
II. n = 1, TO = 1: 



Cc) = 9 \F{iz, k)+g \u{iz, a^, fc), (A.3) 

In the generate cases one obtains: 

a) Ji|i(C,0) = arcsin(^); 

b) Ji|i(C, ^) = ix/3arctan (^^^|^^-arcsin(l-20. 
11. n = 2, m = 1: 



' 4 I 

2 1 ^/-P4(£,£f7) 



te cases one obtains: 



4 I 
g.E{tz,k) 

In the generat( 

a) J2|i(£,0) = l- vT^; 

b) J2|i(£, 
II. n = 3, TO = 1: 



(A.4) 



= -VW^)- f arctan (0i^j . 



Jm{^,^c) = '^-^g-F{iz,k)^ 

' 4 I 

g \Yl{iz, a\ k) - \.J-Pa{^, £c), (A.5) 

In the generate cases one obtains: 

a) J3|i(£, 0) = i arcsin(0 - |£^/^^; 

b) J3|i(£, ^) = ^ arctan (^^|=|) - 

I arcsin(l - 2 * - 1(1 + 20 - 0- 
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II. n = 4, m = 3: 

2{3 + a^)Y 2 1 



^a{9 — cr-') J g I 
C(2(3 + OC + (l~a^)(3-a^)) 
<72(9-a2)V-P4(^,fc) 
In the generate cases one obtains: 



.(A.6) 



a) J4|3(^,0) = 1/^1^-1; 

b) J4|3(C, = ^%/3arctan (^^|^) + 



1 4C^ + 15g + 5 / ^ _ Sttv/S 
18 (?-l/2)^ V 54 " 



Here a := yr^4|^, ^ := a" V^/(l - 0, 

4 



5 := 



V(1 + <t)(3-^)' 



A; := 



l-cr 3 + cr 



l + fj 3-0-' 
and 



1 -cr 



3-C7' 



I ^ ' ' Jo ^{l + x^)il + k^x^) 



-U(iz o? k) ■= f 

i Jo {l+a^x^\/{l+x^){l+k^x-') 



(A.7) 

(A.8a) 

(A.8b) 

(A.8c) 



define for any values < < ^ < 1 tlirec basic elliptic in- 
tegrals in most convenient for us uniform representation. 



Note that for £_c & [O, \/l74] the parameters (A.7) 
have real values: k e [O, l] , .g G [2, 4/\/3] ,ae [O, l/VS] . 

For G ( \/l/4, l) these parameters are complex 
numbers. Nevertheless, for < < ^ < 1 the inte- 
grals (A.2)-(A.6) have real values in this case, too. 

Of course one can use more complicated representa- 
tions of the standard elliptic integrals [7], such that the 
integrals J„|m(Ci^c) will have transparently real values 

forces (^^74,1): 

Taking into account that in this case k = exp{itjj), 
V' € JR, we can perform the transformation of the modulus 

fc ^ fc = i (yk + 1/Vk^ = cos(V'/2) e [0, 1] in all ellip- 
tic integrals. For example, F{z, k) = F{z, k)/2\/k, where 
5-1 = 1 (^zVk + l/zVkj. Unfortunately, such represen- 
tations for E{z, k) and 11(2:, a, k) are quite complicated. 
We will not give them here. They can be obtained com- 
posing several transformations of fc, described in [7]. 
APPENDIX B: TABLE OF THE VALUES OF 
FUNCTION ^r'*(^c) 

TABLE I: The values of function ^r"(^c) 





().()() 


().().") 


0.10 


O.L") 


0.20 


0.2.") 


().:iO 




.9128 






.911.S 


.91()7 


.919:i 


.9.")25 




0.35 


0.40 


0.45 


0.50 


0.55 


0.60 


0.65 




.9564 


.9608 


.9 656 


.9706 


.9756 


.9806 


.9852 




0.70 


0.75 


0.80 


0.85 


0.90 


0.95 


1.00 




.9894 


.9930 


.9959 


.9980 


.9993 


.9999 


1.00 
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